Field induced phase transitions and phase diagrams in BiFeOs — Hke multiferroics 
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The incommensurate magnetic structures and phase diagrams of multiferroics has been explored 
on the basis of accurate micromagnetic analysis taking into account the spin fiexoelecric interaction 
(Lifshitz invariant). The objects of the study are BiFeOz - like single crystals and epitaxial films 
grown on the < 111 > substrates. The main control parameters are the magnetic field, the magnetic 
anisotropy, and the epitaxial strain in the case of films. We predict novel quasi - cycloidal structures 
induced by external magnetic field or by epitaxial strain in the BiFeOz - films. Phase diagrams 
representing the regions of homogeneous magnetic states and incommensurate structures stability 
are constructed for the two essential geometries of magnetic field (magnetic field oriented parallel 
to the principal crystal axis H \\ Cz and perpendicular to this direction HLCz). It is shown that 
the direction of applied magnetic field substantially affects a set of magnetic phases, properties of 
incommensurate structures, character of phase transitions. Novel conical type of cycloidal ordering is 
revealed during the transition from incommensurate cycloidal structure into homogeneous magnetic 
state. Elaborated phase diagrams allow estimate appropriate combination of control parameters 
(magnetic field, magnetic anisotropy, exchange stiffness) required to the destruction of cycloidal 
ordering corresponding to the transition into homogeneous structure. The results show that the 
magnitude of critical magnetic field suppressing cycloid is lowered in multiferroics films comparing 
to single crystals, it can be also lowered by the selection of orientation of magnetic field. Our results 
can be useful for strain engineering of new multiferroic functional materials on demand. 



PACS numbers: 75.85.+t, 75.50.Ee, 75.30.Kz, 75.30.F\. 



I. INTRODUCTION 

Multiferroic materials are compounds that have cou- 
pled magnetic, ferroelectric, and ferroelastic orders. The 
interest in these materials is driven by the prospect to 
control charges by applied magnetic fields and spins by 
applied voltages. They have opportunities for potential 
applications in fields as diverse as nanoelectronics, sen- 
sors, photovoltaics and energy harvesting^"®. 

Although the number of multiferroic materials perma- 
nently increases the one of the most studied multiferroics 
remains the bismuth ferrite [BiFeOj, or BFO). The 
BFO has extraordinary ferroelectric properties^'^'^'^'^, 
a cycloidal magnetic ordering in the bulk^°, and many 
unexpected transport properties such as conductive do- 
main walls^^ or an appropriate bandgap of interest for 
photovoltaics"'^^. It has been used as a blocking layer in 
spin-valves"'^^'"'^^ to control their giant magnetoresistance 
(GMR) by an electric field and also as a gate dielectric 
layer in magnetoelectric field eflFect devices^^. The BFO 
can be interesting for magnonics"'^^ since their magnon 
spectra can be electrically controlled over a wide range^^. 

The BFO has high temperatures of ferroelectric 
Tc — 1083 K and antiferromagnetic ordering 
T/v — 643 K , has an electric polarization of the order 
1 C jm^ and magnetization of the order 5 emujcm^'^''''^^ . 

Since from the sixties the structure and the properties 
of BFO bulk single crystal have been extensively stud- 



jg(j9, 19-21^ Crystal structure comes from the structure of 
ABO3 perovskite oxides. Three types of distortions: rel- 
ative displacement of Bi and Fe ions along < 111 > axis, 
deformations of oxygen octahedral and counterrotation 
of oxygen octahedral around the < 111 > axis reduces 
the perovskite symmetry group to R3c space group. The 
spontaneous polarization caused by the distortions along 
one of eight pseudocubic < 111 > directions. 

Magnetic structure of the BFO in the first approxima- 
tion is G- type antiferromagnet with weak ferromagnetic 
component as it has been established by Kiselev et aP^. 
Further neutron diffraction studies"'^'''^'^'^'' have shown the 
G - type antiferromagnetic structure is subjected to cy- 
cloidal modulation with period of 62 nm lying in the 
plane which contains one of the ferroelectric polarization 
and the propagation vector (Fig. 1). Since eight direc- 
tions of polarization are allowed in the bulk the several 
directions of the cycloid propagation are possible. There 
is a possibility that the cycloid could be either left or 
right handed. However, the cycloids in the bulk were 
found to be of single chirality^". According to^^ the cy- 
cloidal magnetic structure exists below 650 K on cooling 
down to 4 K. The explanation of the complicated spin ar- 
rangement in BFO requires to take into account specific 
spin flexoelectric (or flexomagnetoelectric) interaction. 

The corresponding additional term arising in a free 
energy expansion in crystals belonging to R3c symme- 
try group is known as the Lifshitz invariant. In^^ it 
has been shown that the presence of the Lifshitz invari- 
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FIG. 1. BFO unit cell, schematic illustration of spin cycloid 



ant leads to a solution related to cycloidal spin arrange- 
ment. The relationship between the Lifshitz invariant 
and the Dzyaloshinskii -Morya interaction has been dis- 
cussed in^^. 

The first observations of space modulated structures 
refer to metallic magnets, later on the spiral magnetic 
ordering was discovered in magnetic dielectrics (antiferro- 
magnets), particularly the BFO. At present many mul- 
tiferroics with helical magnetic ordering are known^^'^^. 
The theoretical description of incommensurate super- 
structures in ferromagnetic metals has been elaborated 
by Dzyaloshinskii'^^ . Within the same approach the cy- 
cloidal magnetic ordering in the BFO has been explained 
appealing to the mechanism of inhomogeneous magne- 
toelectric interactions^^"'^^ . The spin cycloid at room 
temperatures is well described by harmonic functions 
sin{kr), cos{kr), where k is the wave vector of spin prop- 
agation. In general case spin distribution in the cycloid 
obeys anharmonical law and is described with Jacobi el- 
liptic functions sn(kr, v), cn{kr, v) where v is the param- 
eter defining the degree of anharmonicity. The change of 
temperature^^, the rare - earth ion doping'^^"'^* , magnetic 
and electric fields^^""*^, stresses induced by orienting sub- 
strate'*^"'*^ are external factors affecting the parameter v. 
The slight structural modifications in the BFO thin film 
can cause drastic changes in the magnetic structure"*^. 
The manifold of magnetic phases are realized in the BFO 
films depending on the type of substrate, the crystallo- 
graphic orientation of the film, the chemical doping, the 
presence of ferroelectric domain structure^^"^^''*^"'*^. 

The aim of presented research is to analyse of possi- 
ble commensurate and incommensurate magnetic struc- 
tures which can be realized in single BFO crystals and in 
the (111) -oriented BFO film, the investigation of tran- 
sitions between modulated and homogeneous magnetic 
states under external conditions. In our consideration 
magnetic field and magnetic anisotropy have been taken 
as the key parameters regulating appearance of magnetic 



phases and their restructuring processes. We reckon to 
the fact that micromagnetic structure is being modified 
under temperature variations of exchange and induced 
anisotropy parameters, and an important factor control- 
ling magnetic states is magnetic field. Our findings show 
that incommensurate magnetic phase is complex, it com- 
prises different phases between which phase transitions 
occur when parameters of a system are changed. Till 
now the incommensurate phase in the BFO - like mul- 
tiferroics has been considered as cycloidal structure with 
spins rotating in the plane passing through the principal 
crystal axis and one of the axes lying in a crystal basal 
plane. We show that new cycloidal phases with three - di- 
mensional spin reorientation arise with change of external 
magnetic field and magnetic anisotropy. The most essen- 
tial geometries of magnetic field (magnetic field oriented 
along principal crystal axis C3 and in the direction per- 
pendicular to C3) have been considered and phase dia- 
grams in terms of magnetic field and magnetic anisotropy 
have been constructed. In frame of the developed model 
one can follow the field and the temperature transforma- 
tions of the BFO micromagnetic structure. 

The paper is organized as follows. In Section II we 
discuss the problem, perform the theoretical model and 
the governing equations. Section III treats homogeneous 
magnetic states, incommensurate states and phase dia- 
grams of the BFO - like multiferroics in the magnetic 
field applied along the principal crystal axis H \\ C3, 
Sections IV, V represent the similar analysis for the 
cases of the magnetic field applied in the basal plane 
of the film. We concentrate on limiting situations con- 
cerning magnetic field oriented along the direction of cy- 
cloid space modulation H \\ OX \\ [110] (Section IV) 
and magnetic field oriented in the perpendicular direction 
H II OY II [112] (Section V), phase diagram related to 
the situation HJ-Cs is presented in Section V. 

Experiments show that cycloid is suppressed in high 
magnetic fields'^^"*^ and is not always observed in the 
BFO thin films*^. The presented in this paper map of 
magnetic states stability (phase diagram) will enable to 
estimate conditions in which definite type of space mod- 
ulated structures exists and also conditions required to 
the space - modulated state destruction. 



II. GENERAL EQUATIONS 

In this section the formalism of micromagnetism ap- 
proach is developed to describe magnetic phases in 
BFO -like multiferroics being considered both in a sin- 
gle crystal and in films. Symmetry and crystallographic 
structure of the film differ from the ones of the single 
crystal and as consequence the crystal and the film can 
possess with different physical properties. The BFO is 
a bright example of the given above assessment. It is 
known that BFO films can demonstrate semiconductor 
properties or even become metallic at definite conditions 
while BFO single crystals are known as insulators. How- 



ever in our consideration we investigate magnetic prop- 
erties of the BFO film limiting a problem with a range 
of parameters doesnt allowing the profound structural 
changes. The (111) -oriented BFO film with rhombohc- 
dral crystallographic structure, the same as BFO crystal 
is considered. 

The determination of magnetic structures and the con- 
struction of corresponding phase diagrams is based on 
the known variational problem of free energy functional 
minimization namely 5^ = 5 J FdV = at the condition 
(5^$ > 0. 

The free energy density of BFO - like crystal is repre- 
sented in a form 



where H is applied magnetic field, the 
anisotropy energy density is represented as 
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where Ku is the constant of uniaxial magnetic anisotropy. 
The effective magnetic anisotropy in films can be dif- 
ferent from the one in single crystals. It is shown below 
that in the case of (111) - oriented BFO films uniaxial 
magnetic anisotropy has additional contribution related 
to the magnctoelastic energy density 
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where 



Fad = AMi • M2 + D[Mx x M2] 



(1) 



(2) 



is the isotropic and the Dzyaloshinskii - Morya exchange 
interactions energy density. Mi and M2 are the sublat- 
tices magnetizations, A is the antifcrromagnetic exchange 
coupling parameter, D = Dric is the Dzyaloshinskii vec- 
tor. Tie is the unit vector oriented along crystal principal 
axis, D is the Dzyaloshinskii parameter. 

In antifcrromagnctism theory it is accepted to use 
ferromagnetic and antifcrromagnetic order parameters 
M — Ml + M2 , L — Afi — M2 or dimcnsionlcss vari- 
ables m = M/2Mo, I = L/2Mo. 

In terms of M, L the Dzyaloshinskii - Morya exchange 
interaction energy density can be rewritten in a form 
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The exchange energy density acquires a form 
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where A is the stiffness constant, li, i = x,y,z are the 
components of the unit antifcrromagnetic vector I, Mq is 
the sublattice magnetization. The non - uniform mag- 
netoelectric interaction energy density known as the Lif- 
shitz invariant is written as 

Fl = Pilx'^xh + ly'^ylz " h'^xlx - Iz'^yly) (5) 

where OX || [110], OY || [112], OZ || [111], the spon- 
taneous electric polarization vector P is supposed to be 
oriented along [111], f3 is the constant of the non - uni- 
form magnetoelectric interaction, its sign is dependent in 
particular on the orientation of vector P. For a definet- 
ness below we suppose that /3 > 0; in the case of the 
BFO multiferroics /3 « 0.6 erg/cm^. 
The Zeeman energy density is given by 
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where B2 is the magnctoelastic constant, uq is the mis- 
match parameter determined over film and substrate lat- 
tice parameters afum, agubs 



Uq = 






(9) 



The lattice mismatch depends on their values in a sin- 
gle crystal, the growth conditions of the heterostructure, 
temperature and thickness of the film'*^-'**. 

As seen from (7), (8) the magnctoelastic interaction 
being taken into account leads to the renormalization of 
the uniaxial magnetic anisotropy constant K^ — K^ + 
B2U0/2. 

At low temperatures T « Tpf {Tn is the Neel temper- 
ature) ferromagnetic and antifcrromagnetic vectors sat- 
isfy to the relations 
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+ m^ = 1 
lm = Q 
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In the relatively weak magnetic fields H << Hex 
{Hex ~ lO^Oe in BFO) these additional conditions al- 
low exclude the vector m from the minimization prob- 
lem, and the free energy density F can be represented in 
terms of the unit vector I and its derivatives. 
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where Heff — MqH + D[l x Bp], Bp — (0,0, 1) is the 
unit vector of spontaneous polarization P oriented along 
the principal crystal axis, h ~ H/Mq. 

Hereinafter we transform to the reduced parameters 
Kc, Km, K.d determined as 
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where Hd = D is the Dzyaloshinskii field, xi- is the 
transversal magnetic susceptibility of antiferromagnet. 



In the calculations carried out below we have chosen 
the values of the parameters characteristic of the mul- 
tiferroics BFO. The literature values of the exchange 
stiffness A is in the range of (2 — 4) • lO^^erg/cm^'^^''*^. 
Below we have taken A — i ■ IQ^'^ erg/cm. The magne- 
tization Mo on various cstimates^^^'^^^^^ is in the range 
(2 — 5)emu/cm^ (films dopped by rare - earth ions have 
larger values of magnetization), we put Mq = 5emu/cm^. 
The field of the Dzyaloshinskii - Moriya interaction, esti- 
mated from measurements of electron paramagnetic res- 



onance^^ is Hd = 1.2-lO^Oe. Transverse susceptibility of 
an antiferromagnet, by ^^ equals x± — 4 • 10~^. Magne- 
tostriction and magnetic anisotropy of the perovskite-like 
antiferromagnets varies quite widely '*^'^^'^^. If we assume 
that the variation of the induced anisotropy in BFO is 
in the range 10^ erg/cm'^ < \Ku\ < 10^ erg/cm'^, then 
it corresponds to a change of the normalized parameter 
Kc in the range —5 < \kc\ < 5. 

The reduced free energy density E — 2AF j (3"^ in terms 
of the variables Kc, 'td, «„, x — xl3/2A acquires the form 
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Due to the identity |Z^| = 1 the vector I may be deter- 
mined by the two coordinates qi = 0,ip (i = 1,2) which 
are the polar and the azimuthal angles in spherical co- 
ordinate frame. The polar angle is measured from the 
equilibrium position of antiferromagnetic vector Iq, the 
azimuthal angle is measured from its projection on the 
orthogonal plane. 

In this paper we restrict ourselves by investigation of 
one - dimensional magnetic structures depending on the 
X coordinate. In this assumption the equation (5$ = 
results in the Euler - Lagrange equation 
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where q'^ = dqi/dx, i ~ 1,2. 

Let qoa{x) be a set of magnetic structures determined 
by equation (13) where a enumerates a set of solutions 
of eq.(13) (a= 1,2,3,...). 

The condition of stability S^^ > of magnetic struc- 
ture qoa arrives to the Sturm - Liouville eigenvalue prob- 
lem 

Ltj{qoa)Sqaj = PaiKSqai (14) 

where the functions 6qai — qi — qoaiix), the differential 
operator Lij and the 'weight' functions pai are fully de- 



termined by the second variational derivative 



Sqi Sqj 



The 



condition (5^$ > requires that all A^ > 0. The corre- 
sponding differential equations will be given below when 
the specific situations are being considered. 

Equation (14) determines a set of the eigenvalues Aq, 
and the eigenfunctions Sqai for the every solution ^oq- 

Parameters Aq, and eigenfunctions Sqai have a simple 
physical sense. The eigenfunctions Sqai can be consid- 
ered as amplitude functions of low energy spin waves or 
magnons. The eigenvalues A^ are proportional to the 
square of magnon frequencies lu related to the qoa ~ raag- 
netic structure namely Aq, = X-L^ai^)/l'^ where k is the 



I 

wave number of magnons, 7 is the gyromagnetic ratio. 
We will use below uj^ik) values instead of Aq. 

It appears to be more convenient to appeal to 
Cartesian representation of antiferromagnetic vector 
I = ilx,ly,lz) for numerical simulation. In this case 
equation (13) acquires a form 



mi 
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where Aq is indetermined Lagrange multiplier. Equation 
(15) can be written as following 
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In its turn equation (14) yields 
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Note that conditions of transitions between magnetic 
phases qa can be described in terms of (14). Eigenvalues 
Aq change with the change of control parameters such as 
magnetic field or intrinsic magnetic anisotropy. In the 
case when the phase qoa loses its stability the parameter 
Aq changes its sign. In other words one of numbers Aq 
becomes equal to zero approaching to the critical point 
of phase transitions. In accordance with Landau theory 
this circumstance determines the soft mode of antiferro- 
magnetic vector oscillations. The condition of the phase 
qoa loss of stability is determined by the vanishing of 
minimal eigenvalues Aq. 

Later on we consider phase diagrams of magnetoelec- 
tric antiferromagnet subjected to magnetic field in terms 
of H and Ke- 
lt should be emphasised that intrinsic magnetic 
anisotropy of BFO - like compounds has a compli- 
cated character, it depends on variations of tempera- 
ture, doping of rare earth ions, stresses arising during 



film growth and lattice mismatch^^'**^"'^^'^^'^^"^'^. The 
consideration and the analysis of physical mechanisms 
giving rise to magnetic anisotropy allow comprehension 
of prevailing factors responsible for magnetic properties. 
Magnetic anisotropy of BFO -like crystals is governed by 
several competing mechanisms including the single -ion 
anisotropy, the anisotropic superexchange coupling, mag- 
neto -dipole interactions. In its turn the single -ion mag- 
netic anisotropy can be divided into several contributions 
attributed to the symmetry of magnetic ions surround- 
ing; the exchange coupling mechanism includes the anti- 
symmetrical Dzyaloshinskii - Morya exchange along with 
other relativistic exchange contributions of quasidipolar 
and non -dipolar types. An overall dependence of mag- 
netic anisotropy of the crystal on internal and external 
parameters, e.g. concentration of rare earth ions and 
temperature variations is determined by the correspond- 
ing behavior of its constituting components^*'^^. 

In respect of BFO -like films an additional contri- 
bution of a surface anisotropy related to a substrate 
should be taken into account. Magnetic anisotropy of 
the film depends on a number of factors: effect of rough- 
ness, the shape of a sample, dipole - dipolar interactions 
g-j-j, 58,59^ Due to the physical origin the surface mag- 
netic anisotropy can be divided in magnetocrystalline 
and magnetoelastic anisotropics^*^. According to the Neel 
model^^ the lack of neighbors at the surface gives rise 
to magnetocrystalline anisotropy so that additional term 
2Ks/t is added to bulk anisotropy where Kg is the surface 
anisotropy, t is the thickness of a film^^'^°"^^. Note that 
such approach fits to ultrathin films where the magni- 
tude of the surface anisotropy is independent on the film 
thickness. Further development of magnetocrystalline 
anisotropy theory in thin films is connected with imple- 
mentation of various methods including first-principles 
calculations®^"^^ , molecular dynamics simulations^^'^^, 
phenomenological models"*^'®^. The other additional con- 
tribution to magnetic anisotropy of films gives the mag- 
netoelastic anisotropy attributed to the magnitostriction 
effect arising due to the lattice mismatch between film 
and substrate. Strains coming from a lattice mismatch 
in epitaxial films distribute from non - magnetic substrate 
into magnetic layer. The coupling between lattice strain 
and magnetization results in magnetostriction inducing 
magnetoelastic anisotropy. In the case of (111) - oriented 
BFO films the magnetoelastic energy density is of a form 
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(18) 



where Bi , B2 are the magnetoelastic coupling coefficients 
accessible from experimental determination^'^'^'' or ab - 
initio calculations^^, l'^ are the components of antiferro- 
magnetic vector, u^^, are the components of deformation 
tensor taken in the Cartesian coordinate frame X' re- 
lated to crystallographic axes [100], [010], [001]. One can 
show that in the coordinate system X connected with 
the principal crystal axis C3 ]|< 111 > the magnetoelas- 



tic energy density given by (18) is rewritten as 

r^ _ B2U0 2 



m.elas. 



(19) 



contributing to the surface magnetic anisotropy. 

This short overview shows that in a frame of a concept 
of magnetic anisotropy one can distinguish films and crys- 
tals by the type of magnetic anisotropy including relevant 
for the considered problem contributions. 

By taking into account the given above consideration it 
is of interest to discuss the change of the ground state of 
the antiferromagnetic multifcrroics by varying the energy 
of the magnetic anisotropy and exchange parameters. We 
treat BFO - like multifcrroics placed in magnetic field 
applied along the principal crystal axis and in the per- 
pendicular direction. In the latter case we investigate 
the influence of variation of the direction of magnetic 
field in the basal plane of a sample relative to the di- 
rection of space modulation of antiferromagnetic cycloid. 
Hereinafter we consider homogeneous magnetic states, 
incommensurate structures and related phase diagrams 
for each orientation of the magnetic field. 



III. MAGNETIC FIELD H \\ OZ APPLIED 
PERPENDICULAR TO THE FILM PLANE 

Let's consider the case H \\ OZ choosing Cartesian 

coordinate frame connected with the principal crystal 

axis C3 II < 111 >: OX || [110], Or || [112], OZ || 
[HI]. 



A. Homogeneous magnetic states 

In the case H \\ OZ the uniform part of the free 
energy density (12) is represented as 



Ea = 7^{Kc + Krnh'')ll = -^{Kc + K,nh^){ll + l^) (20) 



One can see that the magnetic field applied along the 
principal crystal axis renormalizes the constant of mag- 
netocrystalline anisotropy. Dependent on a sign of effec- 
tive anisotropy constant Heff = Kc + K„ih'^ homogeneous 
magnetic state of the "easy plane" |^j,| — 1 or the "easy 
axis" \lz\ = 1 type realizes. As follows from (20) the 
phase \lz\ = 1 exists when h > \/ —nd Hm- The exchange 
and the non -uniform magnetoelectric interactions result 
in the appearance of inhomogeneous magnetic phases and 
the shift of boundaries of homogeneous phase transitions. 
To analyse the stability of homogeneous state it is nec- 
essary to solve eigenvalue problem and to find spectrum 
of natural oscillation frequencies w(A;). For the definite- 
ness we consider non - uniform pertrubation of the ho- 
mogeneous state |/j,| = 1. In accordance with (14) the 
stability condition defining natural frequencies of anti- 
ferromagnetic vector oscillations is determined by 
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(21) 



where k is the wave vector of spin waves, kx is the x - 
projection of vector k indicating the direction of spiral 
propagation. 

Similar consideration can be applied for the analysis of 
the stability of any other "easy plane" state, in particular 
\lx\ = 1- However instability in the last case develops in 
OY direction. 

Critical points of transition from homogeneous mag- 
netic state into modulated structure associated with soft 
mode of oscillations (which is attained at zero values of 
minimum frequency) are determined from the dispersion 
equation (21). Eq. (21) yields the saddle dependence of 
the smallest of natural frequencies oscillation branches 
with two minimums in kx direction. The critical field of 
the transition into homogeneous state is determined by 
the requirement of the vanishing of minimum frequency 
when minimum is attained only at one (e.g. positive) 
value k = kx- Following (21) one can define the critical 
field of the transition from homogeneous magnetic state 
into space modulated structure. 
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(22) 



The curve representing dependence Hc{h) determined by 
(22) is shown on phase diagram (Fig. 4) as line 3. 



B. Incommensurate structures 

To consider the structure and properties of inhomoge- 
neous magnetic states we switch to the polar coordinate 
system with the polar axis aligned along the crystal prin- 
cipal axis and rewrite energy density (12) as 
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Corresponding Euler - Lagrange equations are 



(23) 
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In a general case the system of equations (24) allows a 
set of periodical solutions describing possible space mod- 
ulated structures in multiferroics film. Periodical solu- 
tions differ each from the other by space structure (mag- 
netic configuration), areas of stability dependent on val- 
ues of magnetic anisotropy constant, mismatch parame- 
ter, magnitude and direction of applied magnetic field. 



To consider conceivable periodical structures let us 
start from some simple approximations. In the case 
when magnetic anisotropy and magnetic field are ab- 
sent, eq. (24) has the solution 9o = kxX + kyy = {k ■ r), 
if = arctan(fcj,/fcj:) describing harmonic cycloid. In the 
case when only the uniaxial magnetic anisotropy is taken 
into account the anharmonical solution of (24) described 
by elliptic functions is found as 



sn(- X, ly) 



(25) 



where Keff = Kc + Hmh'^, v is the elliptic modulus < 
V < 1 determined from the minimum of averaged energy 
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where K{v), E{v) are complete elliptic integrals of the 
first and the second kind, (/? is supposed to be constant. 
The other solution differing from (25) by the sign has 
not been considered since it is energetically disadvanta- 
geous. The theoretical analysis of the given above equa- 
tions has been done for ferromagnet and antifcrromagnet 
in works^°''^2_ 

A set of periodical solutions in magnetoelectric an- 
tiferromagnet belonging to the space symmetry group 
RZc = C^y has been considered in^** dependent on the 
constant of uniaxial magnetic anisotropy. It has been 
shown therein that the new type of space - modulated 
structure which is characterized by the conical distri- 
bution of antiferromagnetic vector arises along with the 
plane modulated structure. The first one is denoted here 
as the cone cycloid CC - state which means that spins 
rotate a cone around OY - axis and the second one is 
denoted as the plane cycloid Cy -state which points out 
that ZOX is the plane of spin rotation. Besides that it 
should be noted that the plane modulated structure Cy 
slips into the CC - phase in the case when "hard plane 
anisotropy" attains the critical value corresponding to 
Kcritl = 2.015. 

As was shown at the beginning magnetic field applied 
in [111] direction renormalizcs the constant of magnetic 
anisotropy so spin distribution in the cycloid in this case 
can be also described in a framework of^** by taking into 



account the substitution k 



^eff 



Kc r Kjnll . 



In this item we would like to mention an approach 
allowing estimate the critical field of the transition from 
modulated phase into homogeneous state. It is seen from 
equation (25) that the parameter v defining spin arrange- 
ment in the cyclioid can change in the interval < i^ < 1 
which is accompanied with the subsequent change of the 
effective anisotropy constant < K'^ff — Kc + Kmh'^ < 

Kcrit • 

By taking into account the both relations one can find 
the critical field required for the destruction of space - 
modulated structure following to the condition 



< 



^crit ^c 
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FIG. 3. Spin modulated structure, the scan of projections 
l{x) = {lx{x), ly{x),lz{x)) for the left symmetry CC-solutions, 
solid line corresponds to the dependence lx{x), dashed line 
corresponds to the dependence ly{x), dotted line corresponds 
to the dependence lz{x), Hz ~ l73kOe, Hc = 0.556. 



FIG. 2. Dependence of the period of Cj/-structure on the pa- 
rameter Kc, insert: dependence 6{x), solid curve corresponds 
to Kc = —2.4, dashed curve corresponds to Kc = —2, dotted 
curve corresponds to Kc = —1. 



The plot illustrating dependence of the period of in- 
commensurate structure A = AK{v)v / ^^ K^^r on the ef- 
fective constant of magnetic anisotropy Kc is represented 
in Fig. 2. It is seen that the spiral period changes with the 
varying of effective magnetic anisotropy. The period of 
spiral increases and tends to infinity at the critical value 
of Kc corresponding to Kcrit2 — —2.467. In this case spi- 
ral state disappears, domain walls diverge to infinity and 
the transition into homogeneous "easy axis" state takes 
place. With Kc decreasing the spiral length shrinks, but 
when Kc changes its sign the spiral period increases again 
with the growth of modulus \kc\ and when Kc attains 
values corresponding to Kcriti = 2.015 domains has no 
time to be formed and commensurate structure trans- 
forms into conical state. With further change of effective 
magnetic anisotropy the cone converges to homogeneous 
"easy plane" state at Kcrits = 4'^'*. The scan of I pro- 
jections in CC - phase is shown in Fig. 3. As it is seen 
in Fig. 3 the antiferromagnetic vector goes out from the 
rotational plane in the CC - phase. 



tinguished: homogeneous magnetic states of "easy plane" 
type \ly\ — 1 hereinafter denoted as EP phase and "easy 
axis" type j/zj — 1 denoted as EA phase, two types of in- 
commensurate structures: the plane cycloid Cy and the 
conical cycloid CC being described by the corresponding 
solutions of the Euler - Lagrange equations (24) . The 
Cy - solution corresponds to the cycloid developing in 
ZOX plane, z - and x - components of antiferromag- 
netic vector of Cy phase are described by elliptic Jacobi 
functions, y - component of antiferromagnetic vector 
in this state is equal to zero. The plane of spin rota- 
tion in CC phase is different from ZOX basal plane, 
all the components of antiferromagnetic vector of CC - 
solution are different from zero. Magnetic states contin- 
uously transform each to the other. The plane cycloid 
Cy continuously transforms into the conical cycloid CC 
with the right symmetry of spin rotation, the conical cy- 
cloid CC transforms into the "easy plane" \ly\ =1 state 
when the magnitude of magnetic field and the magnetic 
anisotropy constant enhance. It should be noted here 
that the transition from the modulated Cy state into the 
homogeneous "easy axis" EA phase \lz\ =1 goes over 
the nucleation of domain structure, the transition from 
plane - polarized Cy phase into the conical space - mod- 
ulated structure CC and the transition into homogeneous 
"easy plane state" EP occur along the 2"^^ order transi- 
tion line. 



C. Phase diagram H \\ C3 

The considered analysis together with the computer 
simulation allows reveal a set of magnetic phases realizing 
in multiferroics film in the magnetic field oriented in [111] 
direction coinciding with the principal crystal axis C3. 
The obtained results are presented in terms of the phase 
diagram or the map of incommensurate states stability 
shown in Fig. 4. 

Let us discuss the basic elements of the diagram: the 
possible microstructures, lines and areas of their exis- 
tence and stability. Four types of magnetic states are dis- 



IV. MAGNETIC FIELD H \\ OX APPLIED IN 
THE FILM PLANE 

A. Homogeneous magnetic states 

Let us turn to thhe magnetic field applied in the 
H II OX II [110] direction. We start from determi- 
nation of possible homogeneous magnetic phases which 
can be found out from the uniform part of the free energy 
density 
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Kcl'z "I" n^rnh- 'x + y l^dl^mhiy 



(26) 
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FIG. 4. Phase diagram of a (lll)-oriented BFO film, 
H II C3. Line 1 corresponds to the transition from the 
plane cycloid Cy phase into the "easy axis" EA phase go- 
ing over the phase domains growth, line 2 corresponds to the 
loss of EA phase stability, line 3 corresponds to the 2"'' or- 
der phase transition between the Cy and the cone cycloid CC 
phases, line 4 corresponds to the 2"'' order transition between 
CC and "easy plane" EP phases, the area restricted by lines 
1, 2 is the metastable area of EA and Cy phases coexistence. 



As follows from (26) the symmetrical phase 
I — (0,-1,0) satisfies to minim.um energy condition at 
positive values oi hx > Q at Kc > —^KdKmh. In the case 
^c < —y/i^^di^h the tilted phase Iq — (0, — sin^O: cos^o) 
where sin^g = y/KdKmh/\Kc\ possessing with the energy 
E = KdKmh'^ /"^Kc arises. 

Therefore the transition between symmetrical and 
tilted phases occurs at Kc = —y/Hd^mh in the case 
when non-uniform contributions to the free energy are 
neglected. 

To determine the boundary of the phase transition 
from the symmetrical "easy plane" phase Iq = (0,-1,0) 
into the space - modulated structure we refer to analysis 
of the stability of antiferromagnetic spin structure exist- 
ing in the space uniform state by means of (14) resulting 
in the equation 



UJ - hy/KdKm - K^h - k ) = 4fc^ 



(27) 



allowing determine soft modes of spin excitations indicat- 
ing on a possibility of phase transition. The soft mode 
of the transition corresponds to zero frequency of spin 
oscillations. As in the case considered in Section III A 
one can find that (27) yields the saddle dependence of 
the smallest of natural frequencies oscillation branches 
with two minimums in k^ direction. In the case of soft 
oscillation mode when the frequency tends to zero the 
following condition are to be satisfied 

kx + kl (kc + 2hy/KdKm + Kmh^ - 4) + 

(Kc + h^jKdUm) (Kmh"^ + hy/KdUm) = (28) 

here kx is the x -projection of magnon wave vector fc 



indicating the direction of spiral propagation. The con- 
dition of merging two wave number values corresponds 
to the critical magnetic field value. Therefore the criti- 
cal field of the transition from the homogeneous magnetic 
state into the space modulated structure is defined as the 
minimum positive root of the equation 

(kc -f 2hy/KdKrn + H-mh^ - 4) - 

4 {k.c + hy/KdHjn) {Hmh^ + h^ KdKm) = (29) 

Equation (29) determines the line of the transition 
from the symmetrical phase into the Cy - modulated 
structure corresponding to curve 1 on the phase diagram 
shown in Fig. 10. 

The condition of the transition from the tilted mag- 
netic phase into the incommensurate structure according 
to (14) is of the form 



(w^ - hy/KdUm sin 00 + , 

{oJ^ — h^KdUm sm 9o + Kc cos^ 00 

Akl s\T? 6 



fc^ - 



(30) 



where ^o determines the polar angle of antiferromag- 
netic vector canting in the tilted phase. By tak- 
ing into account that sin^o = —y/KdKmh/nc when 
< h < — Kc/^KdKmh we reduce (30) to 



U! + Kc 



h Kd^r^ 



k^ I X 



(w -\- He — Krnh 



4fc; 



tl KdKr: 



(31) 



Proceeding in a similar way as in the previously con- 
sidered case we find from (31) that the critical field of 
the transition is determined by the formula 

K^ -4 ^c/y/^ 

\/2kI + 8 + 8Kc - KcKd - IQKd/Kc ~ KlJKd 

(32) 



B. Incommensurate structures 

The system of Euler -Lagrange equations determining 
possible incommensurate phases for the field H \\ OX 
is written in a form 

A6' -I- 2 sin^ 6* ( sin .^-^ - cos (^-^ 
\ ox ay 

sin 9 cos 9 [hc — (V(/3)^ — k^/i^ cos^ (/?) — 

y'KdKmh cos 9 sin (fi = 0, (33a) 



V (sin^ 9Wip) 



d9 



d9\ 



2 sin 9 cos (jS^^ — sin (yS-r:^ 
V oy dxj 

Kmh sin sin ip cos if — ^Kd^mh sin 9 cos Lp = (ij33b) 




FIG. 5. Space distribution of the antiferromagnetic vector 
in the cone cycloid CC+, sohd line corresponds to the de- 
pendence lx{x), dashed line corresponds to the dependence 
ly(x), dotted line corresponds to the dependence lz(x), Hx = 
-70kOe, Kc = 0.556, Kd = 0.556. 



accompanied with hysteresis. Such change of the charac- 
ter of the phase transition can be caused as by the change 
in the uniaxial magnetic anisotropy in films and also by 
temperature variations of magnetic parameters in single 
crystals. 

Due to the axial symmetry the space modulation in 
spin subsystem can develop in any direction in the plane 
of the film in the absence of magnetic field. We have con- 
sidered the case when the magnetic field is applied along 
the direction of space modulation of antiferromagnetic 
structure. However the situation when the magnetic field 
is applied at an angle to the direction of space modula- 
tion is possible as well. Below we consider the limiting 
case corresponding to the magnetic field oriented in the 
direction transverse to the direction of space modulation. 
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FIG. 6. The dependence of the space averaged projection of 
the antiferromagnetic vector < /,, > in the cone cycloid CC+ 
on magnetic field starting from the initial plane cycloid Cy 
state, in the CC+-state for the fixed value of Kc = 0.556, Kd = 
0.556, Km = 2.28 • 10~^. Insert: the dependence of the space 
averaged projection of the antiferromagnetic vector < Zj, > on 
the magnetic field calculated for Kc = 2.356. 



In the absence of the magnetic field the only plane cy- 
cloid Cy - state is realized. Numerical analysis of equa- 
tions (33) shows that cone cycloids CC+, CC- differing 
by the direction of spin rotation appear when magnetic 
field is applied. For example the spatial dependences of 
vector I projections in the CC+ ~ state are presented in 
Fig. 5. The both CC+ and CC- states continuously arise 
from the Cy - phase when the magnetic field is applied 
and continuously transform into the homogeneous "easy 
plane" EP state when the magnetic field grows as shown 
in Fig. 6. 

As seen from plots in Fig. 6 and in the insert to this 
figure transitions between CC+ - and CC- -structures 
can be of the f* and the 2"^^ type dependent on the 
value of the reduced anisotropy constant Kc- In the 
area Kc < 2.015 the transition between conical struc- 
tures CC+ and CC- is of nonhysteretic character, in the 
area Hc > 2.015 the transition between conical modu- 
lated structures becomes the first order phase transition 



V. MAGNETIC FIELD H \\ OY APPLIED IN 
THE FILM PLANE 

Consider now magnetic phases and phase transitions 
in the BFO film subjected to the magnetic field applied 
in the H \\ OY \\ [112] direction 



A. Homogeneous magnetic states 

As in previous cases we begin with exploration of ho- 
mogeneous magnetic phases. The uniform part of the 
free energy density (12) acquires a form 



£^0 



1 



-dl 



1 






-\/l^dl^rnlT'>'X 



(34) 



We regard positive values oi hy > 0: in the case Kc > 
the minimum of the free energy density (34) corresponds 
to the symmetrical "easy plane" phase I = (1, 0, 0) (EP), 
in the case Kc < at \kc\ > \/HdKm the tilted phase 
I = (sin6'o, OjCos^o), sin 6*0 = ^/ndKmh/ \kc\ (T) arises. 

The analysis of the stability of the symmetrical phase 
(EP) and the tilted phase (T) are quite similar to the 
case H \\ OX so for this situation we refer to Section 
IV A. The action of the magnetic field applied in OX, 
OY directions in the basal plane due to the axial C3 
symmetry should be equivalent, however the selected di- 
rection of the spiral propagation breaks the symmetry. 
In the case H \\ OX the plane cycloid Cy develops in 
OX direction, in the case H \\ OY the modulated in the 
plane ZOY Cy structure propagates in OY direction. 

By following to such arguments and the calculations 
described above we conclude that the critical magnetic 
field governing the stability of the symmetrical EP phase 
is determined by equation (29) deduced from the equa- 
tion identical to (28), the transition from the T phase into 
the Cy - structure is described by (32) following from the 
equation analogous to (31), in which k^ - projections of 
the vector k are substituted by ky - projections. 
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B. Incommensurate structures 

Turning next to inhomogeneous spin structures arising 
in the magnetic field h — {0,hy,0), we come back to 
the system of Euler - Lagrange equations which in the 
considered case is represented in a form 

A9 + 2 sin^ 9 ( sin 93-^ - cos ip^ 



dx 



2 „• 2 , 



sin^cos^ (kc — (V(p) ~ Kmh sin 

\/HdKmh cos 9 cos if = 0, 



^ 



(35a) 



V (sin^ 9Vip) + 2 sin^ 9 i cos ip—z - sin ip—z ] - 
^ ' \ oy ox J 

Kjnh sin 9 simp cos (fi — y/KjK^h sin 9 simp — 0{35h) 

Periodical solutions of (35) equations describe incommen- 
surate spiral structures differing by magnetic configura- 
tions determined by spatial dependences of 9, ip variables. 
Let us consider solutions with the determined plane of the 
rotation which position is found at the following restric- 
tion 



^=0, 



sin ip ~Q 



(36a) 



(36b) 



Eq. (36) implies spins to be rotated in ZOX plane. The 
law of the spin distribution in the plane cycloid Cy is 
derived from the equation 



d9_ 

dx 



-|- Kc sin 9 + 2y/KdKmh sin 9 — c (37) 



By integrating equation (37) one can obtain 

961 
dx 



±Jc~ KcSin 9 — 2y/KdKmhsin9 (38) 



which can be also represented as 
where 



sm( 



(39) 



72-1 



y C72 - 2^KdKmh-f - Kc 

Y C72 - 2'^^KmKdh - Kc 
C"f = 2y/KdKmH"f^ + 1) + 27KC 



The integration constant c is determined from the min- 
imum condition of the Cy -phase energy 

1 / f)9 \ '^ f)9 1 
£'=-( — 1 - — - -Kc sin"^ 9- y^Kjn^hsin 9 (40) 



averaged over a unit volume by use of (37) 

(E) — ■ / N / a/c— KcSin^^ — 2yKdKrn^sin(^(i0 
^(c) Jo 

27r 1 
— — c 



A(c) 2 



where 



A = 



dO 



KcSin 9 — 2y/KdKmhsin9 



(41) 



(42) 



is the spiral period. 

We exclude the solution (38) with the negative sign 
since it is energetically disadvantageous. Minimization 
of the function (41) with respect to unknown parameter 
c d (E) /dc — leads to the condition 



c— Kc sin 9 — 2^ KdKmh sin 9d9 = 27r (43) 



which allows to calculate the averaged energy of Cy - 
phase (E) = — c/2. Formula (39) together with the con- 
dition (43) allows determining cycloidal structure. The 
scan of cycloid corresponding to Kc ~ 0.556, Kd ~ 0.556 
dX H = 70kOe is represented in Fig. 7. The field depen- 
dence of the cycloid period calculated by use of (42) for 
Kc = 0.556 is shown in Fig. 8. It is seen that there is the 
critical field, the cycloid period increases without limit 
approaching to the critical field value. The unlimited 
growth of the period of cycloidal structure corresponds 
to the transition into space - homogeneous state. 

One can estimate the critical points of the transi- 
tion into homogeneous tilted and symmetrical phases 
'*c — i^c{h) by comparing the energies of correspond- 
ing states. Let us restrict ourselves with negative val- 
ues of uniaxial anisotropy constant since as was shown in 
item A (Section IV) the tilted phase realizes in the area 
|kc| > y/KdKmh^ Kc < 0, d ~ y/KdKmh . As was showu in 
the very beginning in the case d < ~Kc the tilted phase 
possessing with the energy E — d^ /2kc is stable. It yields 
c — —cP' I Kc and the equation determining the critical line 
of transition is represented as follows 



2^ / d^ 

\ Kc sin 9 — 2d sin 9d9 = 27r 



(44) 



In the case d > —Kc the symmetrical phase 9 = -k 12 
possessing with the energy E = —Kc/2 — d occurs, c = 
2d + Kc and the line of the transition from the Cy - 
modulated structure into the symmetrical phase is found 
as 



27r 



^2d (1 - sin 9) + Kc{1- sin^ 9) 



d9 = 27r 



(45) 



However conditions of the transition of the considered 
plane cycloidal structure into homogeneous state are not 
physically meaningful. Homogeneous antiferromagnetic 
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FIG. 7. Space distribution of antiferromagnetic vector in the 
Cy- state, solid line corresponds to the dependence Ixix), 
dashed line corresponds to the dependence ly{x), dotted line 



corresponds to the dependence I, 
0.556. 



■Xx), Hy = 70kOe, Kc = Kd 



state can not be considered as the ground state on lines 
determined by (44), (45) where the period of cycloidal 
structure grows without limit. In the vicinity of these 
lines the ground state is represented by the space cy- 
cloid of conical type modulated along the direction of 
the applied magnetic field. In this connection we made 
numerical analysis of the plane and the cone cycloid ener- 
gies dependent on the variations of magnetic field at the 
different values of reduced anisotropy parameter. Simu- 
lation shows that at low magnetic fields in the restricted 
range of anisotropy energy values the ground state cor- 
responds to the plane cycloid modulated in the direction 
perpendicular to applied magnetic field. With magnetic 
field enhancement the energy of the cone cycloid modu- 
lated in the direction transverse to the direction of mag- 
netic field approaches to the energy of the plane cycloid 
modulated along the magnetic field (Fig. 9). At the crit- 
ical field value the energy of the plane cycloid becomes 
larger than the energy of the cone cycloid modulated in 
the direction perpendicular to the applied magnetic field. 
These values determine the line of the 1"* order phase 
transition which is attained before reaching the critical 
values of magnetic field when the unlimited growth of the 
plane cycloid period occurs. The line of this transition is 
shown by lines 6, 7, 7' on the diagram Fig. 10. 




Hy, kOe 



FIG. 8. Dependence of the period of CC-state on magnetic 
field, Kc — 0.556. The period increases with the growth of 
magnetic field tending to the infinity at the critical field cor- 
responding to the transition into homogeneous "easy plane" 
state. 
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FIG. 9. Dependence of the average energy density of the 
structure on the magnetic field, solid line corresponds to the 
homogeneous _EP-phase, dashed line corresponds to the plane 
cycloid Cy, dotted line corresponds to the cone cycloid CC, 
Kc ~ 0.556. 



C. Phase diagram HJ-Cs 

We summarize the results of the analysis of homoge- 
neous and incommensurate magnetic phases for the con- 
sidered antiferromagnetic system in the BFO film with 
hamiltonian (12) for the case of magnetic field applied 
perpendicular to the principal crystal axis HJ-C3. 

For the dcfiniteness we consider the diagram corre- 
sponding to the magnetic field oriented along OY \\ [112] 
axis, the similar analysis is relevant for the magnetic field 
oriented in the perpendicular direction OX \\ [110]. The 
areas of the existence and the stability of possible ho- 
mogeneous (tilted phases r_|_, T_, "easy plane" phases 
_EP_i_, EP^) and incommensurate magnetic structures 



(plane cycloid Cy, conical cycloids CC+, CC-) arc dis- 
tinguished on the diagram. In the phase T+ ^j, > and in 
the phase T_ ly < 0. The phases EP+ {H > 0) and EP^ 
{H < 0) have the oppositely directed antiferromagnetic 
vectors in the film plane perpendicular to magnetic field. 
The conical phases CC+ and CC-^ differ by the sign of 
antiferromagnetic vector projection on the perpendicu- 
lar to cycloid modulation direction OY. The area of the 
ground states T+ and T_ existence is restricted by lines 1, 
3 and 1', 3'. Lines 3, 3' and 5 are lines of the loss of stabil- 
ity of these ground states. Lines 2 and 3 as well as lines 2' 
and 3' are lines of the 2"'' order phase transition between 
"easy plane" EP+, EP- phases and space - modulated 
structures of the cone type CC+, CC^. They restrict 
areas of the stability of homogeneous symmetrical "easy 
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plane" phases EP+ and EP-. Lines 6, 7, 7' are lines of 
the 1** order phase transition between the plane cycloidal 
space - modulated structure Cy and cone cycloidal struc- 
tures CC+ , CC- . The direction of space modulation is 
perpendicular to the orientation of applied magnetic field 
in the considered Cy structure and parallel to magnetic 
field in CC+, CC- structures. In the area situated inside 
these lines the ground state of niultiferroics corresponds 
to Cy structure, in the outside area the ground state 
corresponds to CC+, CC- structures. Lines 4 and 4' 
are lines of the loss of the stability of the cone cycloidal 
structure with the opposite projection of antifcrromag- 
netic vector on the film plane comparing to the ground 
space - modulated state. They restrict metastable areas 
of the corresponding phases when magnetic field changes 
its sign. In the area restricted by these lines there is 
the line of the 1** order phase transition delimiting CC-^- 
and CC- phases which begins and ends with the tri- 
critical points of the I''* order "liquid - vapour" ~ like 
transition. 

It should be noted here that dashed lines 1, 1' and 7, 
7' are of the qualitative character which is attributed to 
the difficulties of their precise definition in framework of 
numerical analysis. On the lines 1, V unlimited growth 
of the period of space - modulated structure occurs. The 
transition from the cone cycloids CC+ , CC- into the ho- 
mogeneous phases T_|_ , T- goes over the unlimited expan- 
sion of the tilted phase domain. In a vicinity of the f** 
order transition lines 6, 7, 7' the metastable areas of the 
plane cycloid Cy and the cone cycloids CC+, CC- exist, 
however to determine the boundaries of their existence 
the additional analysis is required. 
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FIG. 10. Phase diagram of a (111) -oriented BFO film, 
H J-Cs- Lines 1, 1' correspond to the transition from cone 
cycloids CC+, CC- into tilted phases T+, T-, lines 2, 2' 
restrict the area of homogeneous phases EP+, EP- stabil- 
ity, lines 3, 3' correspond to the 2"'* order transition between 
"easy plane" states EP+, EP- and tilted phases r+, T-, lines 
4, 4' are the lines of the loss of stability of cone cycloids CC- , 
CC+ , line 5 corresponds to the line of the loss of stability of 
T+ and T- phases, lines 6, 7, 7' correspond to the 1"* order 
transition between Cy and CC+, CC- phases. 



VI. CONCLUSION 

Our findings show that the structure and types of in- 
commensurate phases in BFO - like multiferroics, the 
character of phase transitions, phase diagrams substan- 
tially depend on the external magnetic field and the mag- 
netic anisotropy. We stress as well the role of the strain 
induced magnetic anisotropy. A primary sequence of 
phases: homogeneous magnetic state - incommensurate 
phase - domain structure is driven by the constant of 
magnetic anisotropy which is connected with tempera- 
ture, rare earth ion doping, magnetostriction attributed 
to the lattice mismatch between film and substrate. 

It has been shown that cycloidal states in the BFO - 
like multiferroics can be transformed into the trans- 
verse conical spiral structure under the action of uniaxial 
stresses, driving magnetic field or temperature variations. 
Phase diagrams or maps of magnetic phases determining 
the ground state of multiferroics have been constructed 
for the magnetic fields applied along the principle crys- 
tal axis and in the basal crystal plane. These diagrams 
can be used as practical tools to interpret experimental 
data, for strain engineering design in the (111) - oriented 
BFO films with compressive (corresponds to the left part 



of the diagram {kc < 0)) and tensile (corresponds to the 
right part of the diagram(Kc > 0) deformations. 



Another important aspect of the performed research is 
a consideration of the critical magnetic field of the tran- 
sition into homogeneous magnetic state. It is known that 
magnetic field can suppress cycloid but the required de- 
struction value of the field is too high in bulk materials 
which makes them difficult to use. We have shown that 
in epitaxial multiferroics films the destruction field can 
be lowered on account of magnetoelastic energy. In the 
case when magnetoelastic contribution is not sufficient 
the cycloid can be suppressed by the magnetic field which 
critical value is lower than the one in bulk materials. The 
critical magnetic field depends on the direction of the ap- 
plied mgnetic field; the given results show that its value 
becomes lower in the case when the magnetic field is ap- 
plied in the film plane. Our approach gives an oppor- 
tunity to explain experimental observations of antiferro- 
magnetic restructuring of the BFO films not only with 
a magnetic field, but also with an electric field and with 
the change of the temperature in a variable range'*"''*^ . 
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